Elastic Correlations in Nucleosomal DNA Structure 
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The structure of DNA in the nucleosome core particle is studied using an elastic model that 
incorporates anisotropy in the bending energetics and twist-bend coupling. Using the experimentally 
determined structure of nucleosomal DNA [T.J. Richmond and CA. Davey, Nature 423, 145 (2003)], 
it is shown that elastic correlations exist between twist, roll, tilt, and stretching of DNA, as well 
as the distance between phosphate groups. The twist-bend coupling term is shown to be able to 
capture these correlations to a large extent, and a fit to the experimental data yields a new estimate 
of G — 25 nm for the value of the twist-bend coupling constant. 
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The DNA in cukaryotes is tightly bound to an equal 
mass of histone proteins, forming a repeating array of 
DNA-protein complexes called nucleosomes [l[ . A stretch 
of 147 base pair (bp) DNA is wrapped in 1.84 left-handed 
superhelical turns around the histone octamer that forms 
the nucleosome core particle, which is connected via a 
linker DNA to the next core particle. Each nucleosome 
core is a tiny sized spool with a radius of 5 nm and a 
height of 6 nm [2j. The wrapped DNA- histone octamer 
complex is essentially ubiquitous in nature and has a ma- 
jor role in many cell life processes such as gene expression 
and transcription 

In a recent high precision measurement, Richmond and 
Davey have determined the structure of the 147 base pair 
DNA in the nucleosome core particle with 1.9 A resolu- 
tion |4j. They have observed that the structure of the 
bent DNA segment is modulated in the curvature, roll, 
and tilt, and that the twist structure appears to be most 
affected by the specific interactions with the protein sub- 
strate. This experiment provides a wealth of informa- 
tion about the conformational structure of such highly 
bent and strongly interacting DNA, among which we can 
highlight a number of quantitative observations: (1) the 
period of modulation in curvature is set by half of the 
DNA pitch ~ 5 bp, where either the major or the minor 
grooves face the histone octamer, (2) roll appears to have 
the main contribution to the curvature, as it is favored 
over tilt by 1.9 : 1, and (3) the DNA segment is stretched 
by about 1-2 bp as compared to its unbent conformation 




The conformational properties of relatively long DNA 
segments, as well as their elastic response to mechan- 
ical stresses such as pulling forces and torques, have 
been successfully studied using a coarse-grained elastic 



that could take into account 
Other approaches for 
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description s a a 

thermal fluctuations 
studying DNA structure include first-principle computer 
simulations 13|, [14J and phenomenological modelings us- 
ing base-stacking interactions [la, LL6[. In light of the 
recent experimental determination of the nucleosomal 



Kinked bending 



FIG. 1: The conformation of a segment of the nucleosomal 
DNA, between the 13th and the 21st base pairs, is compared 
with that of a simple homogenous elastic rod that bends uni- 
formly. The bending angles between two consecutive base 
pairs give a measure of the non-uniformity in the shape of 
the highly bent DNA segment. Note the sharp bends (kinks) 
near the 15th and the 20th base pairs. A schematic definition 
of the bending components Roll and Tilt is presented in the 
right panel. 



DNA structure, a corresponding theoretical analysis is 
called for, and one naturally wonders which of the above 
approaches could more easily accommodate the addi- 
tional complications due to the high degree of bending 
and the specific DNA-protein interactions. 

Here, we attempt to use an augmented elastic descrip- 
tion to account for a number of observations made by 
Richmond and Davey. We consider an elastic energy ex- 
pression that includes anisotropic bending rigidities and 
twist-bend coupling. We show that the anisotropic bend- 
ing elasticity is responsible for the modulations in curva- 
ture with the period of 5 bp [13] , and calculate the shape 
of DNA, an example of which is shown in Fig. [1] Using 
an analysis of the experimental data of Ref. [J] , we show 
that the specific features in twist and bend are correlated 
to a large extent through the twist-bend coupling, and ex- 
tract an estimate of G = 25 nm for the twist-bend cou- 
pling constant that best describes this correlation. We 
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FIG. 2: (a) Roll, (b) tilt, and (c) the super-helix-component of 
the difference of the lengths of the two phosphodiester chains. 
The filled circles are the experimental data taken from Ref. 
[3, the hollow triangles (squares) are calculated points with 
(without) the twist-bend coupling. 



calculate the components of curvature as well as the ax- 
ial strain using the experimental values for twist as input, 
and find that roll is favored over tilt by 1.7 : 1, and that 
the overall stretching of DNA is about 1 bp, both in en- 
couraging quantitative agreement with the experiment. 
We have also studied another parameter called APP a , 
which shows the super-helix-component of the difference 
of the lengths of the two phosphodiester chains and 
found good agreement with the experimental data. 

To study the structure of the nucleosomal DNA, we 
consider a simple model in which the molecule is repre- 
sented as an elastic rod [17| . The rod is parameterized by 



the arclength s and at each point, an orthonormal basis 
is defined with the unit vectors ei(s), e2(s), and e 3 (s), 
where e\ shows the direction from minor groove to ma- 
jor groove, and £3(5) is the unit tangent to the axis (see 
Fig. [T|). Note that due to the helical structure of DNA, 
ei(s) and &i{s) rotate with the helix. The deformation of 
the double helix is characterized by the angular strains 
^1,2 (s) corresponding to bending in the plain perpen- 
dicular to ei^ (s) and ^(s) corresponding to twist and 
torsion. We should mention that each slice of the rod is 
labeled by s, which is corresponding to arc length along 
the unstressed helix axis and so always changes from 
to L. The actual arc length along the deformed axis is 
given by ds' , which in terms of axial strain a(s) one have 
ds' = [1 + a(s))ds [HI]. The elastic energy for the defor- 
mation of DNA in units of thermal energy is written as 
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where A\ and A2 are the bending rigidities for the "hard" 
and "easy" axes of DNA cross section, C is the twist 
rigidity, B is the stretch modulus, D is the twist-stretch 
coupling, and G = Gljq is the twist-bend coupling. The 
spontaneous twist of the helix is defined via its pitch P 
as ujq = 2ir/P, and for B-DNA we have loq = 1.85 nm -1 . 
Note that the twist-bend coupling G can be ruled out 
by symmetry (^ to —tp) for a non-chiral rod u>q = 0, and 
thus its presence is a direct consequence of a spontaneous 
twist structure, as manifested by the form G = Gujq. 

The local strains can be written in terms of the cur- 
vature k(s), the torsion r(s), and the twist angle ip(s) as 
ill = Ksini/;, ^2 = kcos?/>, and ^3 = r + d s ip (jj. It can 
be shown that the mean curvature and torsion imposed 
by the nucleosomal structure are n av = R/ (R 2 + v 2 ) 
and T av = — v j (El 2 + v 2 ), where R is the radius of the 
nucleosome and 2ttu denotes the pitch of the wrapped 
DNA around the histone octamer. For the nucleosome 
case, R ~ 41.9 A and 2niy ~ 25.9 A Q], which yields 
T av/ K av = 0.098. Since typical values of d s ip is of the or- 
der of Wo, we can estimate the relevance of torsion with 
respect to twist by the ratio r av /dj = 0.012. These es- 
timates justify neglecting torsion with respect to twist 
and curvature. After defining A = \{A\ + A2), and 
A' = ^(Ai — A 2 ), we can write Eq. as 
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+ 2G(d s ip — luq)k cosip] . (2) 

For the coupling parameters, we use A = 50 nm [lfj| . 
A' = 30 nm 0, B = 78 nm C = 75 nm [2l|, and 
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D = 15 nm [22|. Since the only estimate available for 
the value of the twist-bend coupling constant G has been 
rather indirect 0] , we treat G as a tuning parameter and 
find its value by fitting to the experimental data of Ref. 

1- 

To find the shape of the DNA, we need to incorpo- 
rate the interactions with the binding substrate. The 
wrapping of DNA by an overall angle of = 2tt x 1.84 
around the protein octamer can be imposed as a global 
constraint, which reads ds(l + a)K — 0. To take ac- 
count of the specific and local interactions, we assume 
that it is mostly the twist degree of freedom that is af- 
fected by those interactions [23j ; a view that is supported 
by the experimental observations of Ref. [1] . This allows 
us to write the interaction energy term, which is the final 
addition to the energy expression for the DNA conforma- 
tion Eq. ([2|), as VdnA— Histoae[^]- We can now minimize 
the sum of -Edna + Vdna— Histonc [V 1 ] with respect to k, a, 
and ip, subject to the wrapping constraint above. This 
yields expressions for the curvature and the axial strain 
as functions of the twist angle as 

_ (d s tp - ljq) [fiu! D + BGujI cost/)] - hBloI 
K ~ [p- — BlOq (A - A' cos 2%p) ' ( ' 

(d s xp — ujq) {Dljo (A — A' cos ip) + fiG cos tp] — 
a = y?-Bul (A - A' cos 2ip) ' 

(4) 

where /i is the Lagrange multiplier for the constraint, 
as well as an equation for the twist that contains un- 
known interaction terms arising from VnNA-Histone [ip] ■ 
Instead of elaborating on the possible forms of this in- 
teraction, we simplify the procedure by directly read- 
ing off the twist angle from the experimental results of 
Richmond and Davey. This allows us to make a direct 
comparison between the calculated values for the cur- 
vature and stretching and the corresponding experimen- 
tally measured values for them, and hence put the elastic 
model to a stringent test. 

The twist angle is obtained from the twist strain by in- 
tegration, and the corresponding integration constant is 
chosen by noting that at the dyad axes of the nucleosome 
the major groove-minor groove direction is perpendicu- 
lar to the surface of the nucleosome core particle. This 
means that the first base pair in either the left or the 
right half of the nucleosomal DNA should have an offset 
twist angle of ipo — 2tt x ^ £ = 36°. Finally, we note 
that following Ref. [3| the calculations are performed for 
half of the DNA length corresponding to 73 base pairs, 
to make a direct comparison possible. 

In Figs. [2]a, and [2j), we have plotted the calculated 
roll K = ±^M7 2 = ^nbcosib and tilt T = i^&fii = 
^■nbsuiip [24| as functions of the position of the base 
pairs, together with the experimental data of Ref. [3], 
where b = 3.4 A is the base-pair step for B-DNA. The two 
quantities appear to be modulated with a near periodicity 




FIG. 3: The Fourier transform n q of the curvature. The filled 
circles are the experimental data taken from Ref. Q, the 
hollow triangles (squares) are calculated points with (without) 
the twist-bend coupling. The peak at q — 15.6 corresponds 
to a periodicity of 5 bp. (A large peak at zero corresponding 
to the average of the curvature is eliminated to enhance the 
resolution of the figure.) 

of 10 bp, which is imposed by the near periodicity in ip. 
The best fit to the experimental data yields G = 25, 
which results in the roll being favored over tilt by 1.7 : 1, 
to be compared with the experimental value of 1.9 : 1. 
This shows that in such a highly bent structure, the DNA 
prefers by a ratio of nearly 2 to 1 to use the bending over 
the easy axis as opposed to the hard one. 

To make a more refined quantitative comparison with 
the experiment, we consider the net curvature (as op- 
posed to its components) and take its Fourier transform 
defined as n q = -^E™=1 K^e 2 *^ 8 - 1 )^ 1 )/" for a list k s 
of length n (n — 73 here), to better resolve its features. A 
plot of the absolute value of the curvature Fourier trans- 
form is shown in Fig. [3l comparing the experimental 
data with the calculated ones. Note that the absolute 
value of the Fourier transform is symmetric with respect 
to the transformation q — > n — q, hence only the first 
half of the plot is shown. The Fourier transform of the 
curvature shows a distinct peak at q = -y + 1 = 15.6 
corresponding to a periodicity of 5 bp, which is a result 
of the anisotropic bending elasticity [l7j. Moreover, Fig. 
|3] shows that while the simple elastic description fails to 
account for the detailed features of the curvature without 
a twist-bend coupling, once equipped with such a term it 
can give a considerably improved account, with the best 
fit corresponding to G = 25 nm. We also note that using 
the values for the curvature the shape of the bent DNA 
can be determined upon integration. We have provided 
such an example in Fig. [TJ where the presence of two 
kinks at the distance of 5 base pairs from each other can 
be visibly noted. 

Using Eq. (j4|), one can also determine the stretch- 
ing of the nucleosomal DNA. In Fig. [H the axial strain 
a is plotted as a function of the base-pair position. A 
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a full microscopic computer simulation of such a large 
DNA-protein complex appears to be out of reach with 
the computational power at hand, such simplified phe- 
nomenological approaches could be helpful in under- 
standing the structural properties of biomolecules. 

We thank R. Bruinsma, H. Flyvbjerg, T.B. Liverpool, 
P. D. Olmsted, Z.-C. Ou-Yang, W.C.K. Poon, M. Rao, 
H. Schiessel, and A. Travers for very helpful discussions. 
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FIG. 4: The calculated axial strain as a function base-pair 
position. While detailed experimental data is lacking for the 
stretching at each base pair, the overall stretching obtained 
by integration over the above data gives 1 bp in agreement 
with the observations of Ref. 



positive (negative) value of a shows a stretching (com- 
pression) for the corresponding base-pair. The overall 
length of the DNA in the nucleosome can be found as 
J L ds (1 +a) = 148 bp, which suggests that the DNA in 
the nucleosome is stretched by about 1 bp, in agreement 
with the observations of Richmond and Davey [4J . 

We have also calculated the difference between the 
components of the phosphate-phosphate distances lying 
parallel to the path of the super helix APP a = (£i — £2) ■ 
es(s), where li gives the phosphate-phosphate distances 
on the zth strand. Using the geometrical definitions, 
we find APP a = cos^ (d s ip) (1 - cos 6ft), where 
(d s ip) 2 , and d = 2 nm is the diameter of the 
The calculated values of APP a are 



n = 

undeformed DNA. 
shown in Fig. [2t, which are in good agreement with the 
experimental values taken from Ref. [4|. 

We have also examined the effect of other elastic 
couplings — such as the bend-stretch coupling — by trying 
to fit to the experimental data, and have found no signif- 
icant effect. Higher order elastic terms such as the cubic 
terms in curvature etc. 0], are expected to add cor- 
rections of the order of nd, which become (marginally) 
important in the highly bent (kink) regions. While the 
addition of such terms would definitely help improve the 
results quantitatively, the fact that it will introduce more 
unknown coupling constants make such a direction not 
particularly appealing. Moreover, there are also other 
structural properties of the nucleosomal DNA such as 
shift and slide 24J, which appear to be beyond such a 
simple elastic description. This suggests that a more 
promising direction for an improved theory that can bet- 
ter describe the conformation of nucleosomal DNA is a 
generalization of the base-stacking model of O'Hern et 
al. ~ 
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In conclusion, we have shown that an elastic theory 
that takes into account anisotropic bending and twist- 
bend coupling can account to a considerable degree for 
the observed structure of the nucleosomal DNA. Since 
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